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We summarize recent higher-order QCD results based, directly or indirectly, on the Mellin-space computation 
of the next-to-next-to-leading-order splitting functions governing the evolution of hadronic parton distributions. 
Specifically, we briefly present third-order results for the coefficient functions in inclusive deep-inelastic scattering 
(including the structure function F3 not published so far), the on-shell quark and gluon form factors, and the 
total cross section for Higgs production at hadron colliders. 



1. Introduction 

For the next decade, the highest-energy ex- 
periments in particle physics will be performed 
at the (anti-)proton-proton colliders Tevatron 
and LHC. The cross section for, for example, the 
inclusive production of a high-p T hadron h at 
such machines can be schematically written as 



'VP' 



■ h+x = f P * /p * off * Dp, 

fJ'J" 



(1) 



Here f p 



stands for the universal momentum dis 
tributions of the partons / in the proton p, f = 
QiiQii 9 with i = 1, . . . , rij, where is the num 
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ber of effectively massless flavours, and Dp, are 



the corresponding fragmentation functions of the 
final-state hadron h. a represents the partonic 
cross section for the process under consideration, 
and power corrections have been disregarded. 

Hence all quantitative collider studies of the 
standard model, and of expected and unexpected 
new particles, require a precise understanding of 
the partonic luminosities and of the QCD correc- 
tions to the hard cross sections. For many impor- 
tant processes, like Higgs-boson production, at 
least the second-order (NNLO) QCD corrections 
need to be included, i.e., the third term in 



a = a™ 
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The consistent inclusion of requires parton 
distributions, and in the special case of Eq. (1) 
also fragmentation functions, evolved with the 
corresponding NNLO splitting functions 

(3) 
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The computation of the third-order splitting 
functions for the (unpolarizcd) parton distribu- 
tions has been completed two years ago [1,2], see 
also Refs. [3,4] . As briefly recalled below, that cal- 
culation was set up such that it could be extended 
to the third-order coefficient functions in inclu- 
sive deep-inelastic scattering (DIS) [5,6]. Here 
we briefly discuss these coefficient functions, as 
well as some further results based on their calcu- 
lation [7-9], see also Refs. [10-12]. A discussion 
of related NNLO results [13] for the evolution of 
fragmentation functions can be found in Rcf. [14]. 

2. Third-order DIS coefficient functions 

A calculation of inclusive DIS, {7, W}*(q) + 
f(p) -» X with Q 2 = -q 2 > and p 2 = 0, 
has been performed in Refs. [1,2,6] up to the 
third order in the strong coupling a s = a s /(4ir). 
The results have been obtained first for all 
even or odd values of the Mellin variable N, 
via the optical theorem and the three-loop for- 
ward Compton amplitudes {7, W}*(q) + f(p) — > 
{7, W}*(q) + f(p). From these expression the 
complete Bjorkcn-x space results have been ob- 
tained by an automated Mellin inversion [15,16]. 
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This approach had the crucial advantage of al- 
lowing to check the extensive and involved FORM 
[17,18] codes, at almost any stage, by falling back 
to the Mincer program [19,20] employed in the 
previous fixed- N calculations of Rcfs. [21-23]. An 
independent check of the 7*-exchange non-singlet 
results for N = 16 has been performed in Ref. [24] . 

The pole terms of the unfactorized coefficient 
functions, supplemented by a corresponding cal- 
culation of DIS by exchange of a scalar </> directly 
coupling to gluons (like the Higgs boson in the 
heavy top-quark limit), deliver the complete set 
of NNLO splitting functions [1,2]. The number 
of relevant diagrams is shown in Table 1. The fi- 
nite pieces of the three-loop 7*/ amplitudes, when 
calculated also for the p^p w projection of the 
hadronic tensor, lead to the a 3 coefficient func- 
tions for Fl and F2 in electromagnetic DIS [5,6]. 

The former coefficient functions complete the 
NNLO description of photon-exchange DIS in 
massless perturbative QCD, see, e.g., Ref. [25]. 
At large Bjorken-x the latter quantities dominate 
the N 3 LO corrections for F2 [26,27], thus facilitat- 
ing improved determinations of a s from data on 
dccp-inclastic scattering. For recent fit-analyses 
see Ref. [28] and references therein. 



The small- a; behaviour of the coefficient func- 
tions for F2 is illustrated in Figs. 1 and 2. The for- 
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Table 1 

Numbers of diagrams employed in the calculation 
of the three-loop splitting functions and the cor- 
responding coefficient functions, h denotes the 
standard ghost used to simplify the gluon polar- 
ization sum. The overall sums include the number 
of Lorentz projections (last column) required for 
extracting the coefficient functions for F% and Fl ■ 



mer shows the new N 3 LO contributions 



x), 



as in Eq. (2) using the expansion parameter a s = 
a s /(47r), together with successive approximations 
by the dominant contributions for x — ► 0. Even 
disregarding the additional effect of the ubiqui- 
tous Mellin convolution, cf., e.g., Ref. [2], the 
leading logarithms do not provide a useful ap- 
proximation in the x-range accessible to collid- 

(3) 

ers. In particular, a 20% accuracy for c 2ns (x) is 
reached with one and two small- a; logarithms only 
at x < 10~ 50 (sic) and x < 10~ 14 , respectively 
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Figure 1. The small- a; behaviour of the third- 
order flavour non-singlet (top) and singlet (bot- 
tom) coefficient functions for the structure func- 
tion F2 in electromagnetic DIS. Also shown are 
the respective approximations obtained by suc- 
cessively including the terms leading for x — > 0. 
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The perturbative stability at small x is illus- 
trated in Fig. 2 for the gluonic contribution at 
Q 2 w 30...50 GeV 2 . The third-order correc- 
tion is far smaller than the second-order contri- 
bution first calculated in Ref. [29], exceeding 1% 
of the leading-order result only at x < 2 • 10~ 5 . 
It should be noted, on the other hand, that the 
perturbative expansion appears to become unsta- 
ble at x < 10~ 4 at very low scales, Q 2 « 2 GeV 2 
[5,6], due to the larger a s and especially the flat- 
ter small- a; shape of the parton distributions. 
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Figure 2. The perturbative expansion up to three 
loops (N 3 LO) of the gluon contribution to the 
flavour-singlet structure function F^ for a s = 0.2, 
n f = 4 and xg = 1.6 x " 3 (1 - x) 4 - 5 (1 - 0.6 x a3 ). 
The results have been normalized to the leading- 
order result given by the singlet quark distribu- 
tion xq s = 0.6x~ a3 (l - x) 3 - 5 (1 + 5.0 x - 8 ). 

The convergence of the perturbation series at 
large- a; / large- TV can be conveniently illustrated 
by displaying the value cta n \N) for which the 
effect of the n-th order term Ca(N) is half as 
large as that of Ca [N). This quantity would 
be order- independent for a geometric scries, while 
a systematic decrease with increasing n would in- 
dicate the asymptotic character of the expansion. 
Si^ris(iV) are shown in Fig. 3 for a = 2, L (note 
the different scales). No sign is observed of an im- 
minent breakdown of the perturbative expansion 
at phenomenologically relevant values of N. 




Figure 3. The TV-dependent values of a s at 
which the effect of the n-th order (N"LO for F2, 
N Il_1 LO for Fl) non-singlet coefficient functions 
is half as large as that of the previous order. 



The third-order coefficient function has also 
been computed for the (odd-iV) non-singlet quan- 
tity F^ +9 , the qW entry in Table 1. The lengthy 
exact results and a full discussion will be pre- 
sented elsewhere 1 [30]. Here- we only provide a 
compact x-space parametrization, which is suffi- 
ciently accurate for all numerical purposes and 
which can be easily transformed to Mellin space 
for use which complex- N codes like Ref. [31] . This 
parametrization is given in Eq. (4) using the ab- 
breviations xi = 1 — x, Lq = In x, L\ = In xi, 
and T>k = [x~^ 1 L , {\ + for the usual -(--distributions. 
The factor /Z02 (= 1 for the numerical evalua- 
tion) indicates the d ahc d a b c contribution entering 
at this order for the first time, cf. Ref. [23]. 
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,(3) 



512/27 V 5 - 5440/27 X> 4 + 501.099 V 3 + 1171.54 V 2 - 7328.45 V x + 4442.76 T> 

- 9172.68 5{ Xl ) - 512/27 L\ + 8896/27 h\ - 1396 + 3990 L\ + 14363 L l 

- 1853 - 5709 x + xx! (5600 - 1432 x) - L L X (4007 + 1312 L ) - 0.463 



293.3 L n 



1488 



496.95 Lg 



4036/81 L\ - 536/405 L 



n f { 640/81 V A - 6592/81 V 3 + 220.573 £> 2 + 294.906 V 1 - 729.359 V + 2575.46 c5(xi) 

- 640/81 h\ + 32576/243 L\ - 660.7 h\ + 959.1 L x + 516.1 + «j (635.3 + 310.4 x) 

- 465.2 x + 31.95 xiL? + L L X (1496 + 270.1 L - 1191 £i) - 1-200 xL% + 366.9 L 
+ 305.32 Ll + 48512/729 L\ + 304/81 } 

n| { 64/81 X> 3 ~ 464/81 V 2 + 7.67505 Pi + 1.00830 V Q - 103.2602 5( Xl ) - 64/81 h\ 
+ 992/81 h\ -49.65 Li + 11.32 + 51.94a; - xx x (44.52 + 11.05 x) + 0.0647 xL% 

- LqL x (39.99 + 5.103 L - 16.30 L{) - 16.00 L - 2848/243 L% - 368/243 L% } 
fl 02 n f {2.147 L\ - 24.57 ii +48.79 - x x (242.4- 150.7 x) - L Li (81.70 + 9.412 Li) 

+ xL (218.1 + 82.27 L 2 a ) - 477.0 L - H3.4 L 2 + 17.26 - 16/27 L 5 } sci . 



Likewise the (experimentally irrelevant) coeffi- 
cient functions for DIS via Higgs exchange have 
been computed to three loops in the heavy top- 
quark limit, see Table 1. The x — » 1 limit of these 
(also still unpublished) expressions forms the ba- 
sis for the Hgg results in the next two sections. 

3. Massless quark and gluon form factors 

The form factors of quarks and gluons are 
gauge invariant (but infrared divergent) parts 
of the perturbative corrections to inclusive hard 
scattering processes. They summarize the QCD 
corrections to the qqX and ggX vertices with a 
colour-neutral particle X of either space-like or 
time-like momentum q. These quantities are also 
key ingredients in the infrared factorization of 
general higher-order amplitudes [32,33]. 

The relevant amplitude for the space-like j*qq 
case, directly entering the photon-exchange coef- 
ficient functions discussed above, reads 



(4) 



The corresponding Hgg vertex defining J- g is an 
effective interaction in the heavy top-quark limit, 



£cff 



1 



f * tt s~i a si a,, 



(6) 



where G^ v denotes the gluon field strength ten- 
sor, and the coefficient Ch includes all QCD 
corrections to the top-quark loop, known to 
N 3 LO [34], see also Refs. [35,36]. 

The quark and gluon form factors were directly 
calculated at two loops in Refs. [37] and [38] to or- 
der e°, respectively, and extended to (all) higher 
powers of e in Refs. [7,8,39]. For the status of a 
direct three- loop calculation sec Rcf. [40] . 



The £" 



e 1 pole terms of the three-loop 



r„ 



(5) 



form factors can be extracted from the third- 
order coefficient functions for DIS [7,8]. For this 
purpose we consider the bare (unrenormalizcd 
and unfactorized) partonic structure functions F b 
for 7*q — > qX and <f>*g — » gX in the limit 
x — > 1. Keeping, at each order a™, only the 
singular pieces proportional to 6(1 — x) and the 
+-distributions 



where e q represents the quark charge and Q 2 = 
—q 2 the virtuality of the photon. The gauge- 
invariant scalar function J- q is the space-like 
quark form factor which can be calculated order 
by order in the strong coupling in dimensional 
regularization with D = 4 — 2e. 



V, = 



ln'(l -x) 
(1-x) 



I = 1, 



2n-l, (7) 



these results are compared to the general struc- 
ture of the n-th order contribution F„ in terms of 
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the /-loop form factors Ti and the corresponding 
real-emission parts Si , 

F b = 6(1 -x) 

= 2Fi8{l-x)+Si 
F 2 b = (2T 2 + T x 2 ) 5(1 - x) + 2F 1 S 1 + S 2 
F 3 b = (2JF 3 + ITxTi) 5(1 - x) 

+ (2JF 2 + :F 1 2 ) 5i + 2^i5a + 5 3 . (8) 

The x-dependence of the real emission factors 
<S/c in inclusive deep-inelastic scattering is of the 
form Sk(fk,e), where the D-dimensional +- dis- 
tributions fk. £ are defined by 

f k , E (x) = + 

= 4*(l-a:)+E^^«- 0) 

i=0 

Thus, exploiting this particular analytical de- 
pendence, the n-loop form factor T n can simply 
be extracted by the substitution 

Vo - — 5(1 - a:) - V , (10) 

2—1 

once the combinations of lower-order quantities 
in Eq. (8) - determined before to a sufficiently 
high order in e - have been subtracted from the 
calculated results for F b . As 5(1 — x) enters with 
a factor 1/e, this extraction loses one power in e. 
Hence from the third-order calculation to order 
e°, as performed for the coefficient functions, only 
the pole terms of J-3 can be obtained in this man- 
ner. So far the e°-term has only been derived for 
the fermionic (n^) part of the quark form factor 
[8] , via extending the corresponding DIS calcula- 
tion [41] by one order in e. 

The reader is referred to Refs. [7,8,11] for a 
further discussion of these results, including their 
very interesting structure in the context of the 
exponentiation of the form factors [42-45]. We 
would like to note, however, that the a 3 e _1 coef- 
ficients of the highest C-function weights, C2C3 an( i 
£5, agree with the results inferred, using the con- 
jecture of Ref. [46], from the calculation in Af = 4 
Super- Yang-Mills theory in Ref. [47]. 



4. Higgs productions at (almost) N 3 LO 

Armed with the third-order splitting functions 
[1,2] and the results of the previous section, one 
can now derive all soft-enhanced (-(--distribution) 
contributions to the N 3 LO cross sections for 
lepton-pair and Higgs boson production at col- 
liders [9,50], see also Refs. [51,52]. Analogous to 
Eq. (8), the soft limit of the bare cross sections 
W h for qq — > 7* — > l + l~ and gg — ► H reads 

Wo = 6(1 -x) 
W\ = 2 Re T x 6(1 - x) +5i 
W% = (2RcT 2 + \T 1 \ 2 )5(l-x) + 2RcT 1 S 1 
+ S 2 

= (2 Re T 3 + 2\T X T 2 |) 6(1 - x) (11) 

+ (2 Re T 2 + \T x f)S x + 2 Re^i5 2 + S3 , 

where, of course, T now denotes the time-like 
quark or gluon form factor, known by analytic 
continuation from q 2 = — Q 2 < to q 2 > 0. 

The real-emission contributions Sk depend on 
the scaling variable x = M 2 , H /s. Here the de- 
pendence of Sk on x is of the form Sk(f 2 k.e), i-e., 
00 

Sk = hk,e 2fcs M e ' • ( 12 ) 

l=- 2 k 

With the known time-like form factors, the ex- 
pansion coefficients s k l of the soft function S k can 
be derived recursively as far as they are subject to 
the KLN cancellations and the mass-factorization 
structure relating the remaining poles to the 
splitting functions (3). Employing the results 
of Refs. [7,8] and [1,2], the third-order terms 
s 3 _ 6 . . . 53,-1 can be obtained. Due to Eq. (9) 
this is sufficient to derive all -(—distribution con- 
tributions to the third-order coefficient functions. 
The explicit results can be found in Ref. [9]. 

Focusing on Higgs production, we note that a 
good approximation (to about 10% or less) to the 
double convolutions g * g * a is obtained at NLO 
and NNLO [53-55] by transforming to A-spacc 
and keeping only the ln fe A and A terms arising 
from the -(--distributions in and &( 2 \ Con- 
sequently the above results facilitate a sufficient 
approximation to the complete N 3 LO correction, 
with a (conservative) error estimate of 20%. 
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The resulting predictions arc illustrated for the 
LHC in Fig. 4, where all higher-order contribu- 
tions have been calculated in the heavy top-quark 
approximation, but are normalized to the full 
lowest-order cross sections. Considering these 
and other results [9], 5% at the LHC, and 7% 
at the Tevatron, appears to represent a con- 
servative estimate of the improved cross-section 
uncertainty due to the truncation of the pertur- 
bation series at the (approximated) N 3 LO. 
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Figure 4. Pcrturbative expansion of the total 
cross section of Higgs production at the LHC for 
the parton densities of Refs. [48,49]. Top: depen- 
dence on the Higgs mass Mh- Bottom: renor- 
malization scale dependence for Mh = 120 GeV. 



5. Summary 

We have computed the third-order coefficient 
functions for the most important structure func- 
tions in deep-inelastic scattering [5,6,30]. This 
first calculation of the three-loop corrections to 
one-scale partonic cross sections offers new in- 
sights into the process at hand, and facilitates 
improved determinations of the strong coupling 
constant a s . Moreover it has lead to further im- 
portant third-order (and all-order resummed, cf. 
Ref. [10,11]) results for the on-shell quark and 
gluon form factors [7,8] and the cross section for 
Higgs boson production at proton colliders [9] . 
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